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RADIO IMAGING OF THE NGC 1333 IRAS 4B REGION 
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ABSTRACT 

The NGC 1333 IRAS 4B region was observed in the 6.9 mm and 1.3 cm continuum with an angular res- 
olution of about 0.4 arcseconds. IRAS 4BI was detected in both bands, and BII was detected in the 6.9 mm 
continuum only. The 1.3 cm source of BI seems to be a disk-like flattened structure with a size of about 50 
AU. IRAS 4BI does not show any sign of multiplicity. Examinations of archival infrared images show that the 
dominating emission feature in this region is a bright peak in the southern outflow driven by BI, corresponding 
to the molecular hydrogen emission source HL 9a. Both BI and BII are undetectable in the mid-IR bands. The 
upper limit on the far-IR flux of IRAS 4BII suggests that it may be a very low luminosity young stellar object. 
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1. INTRODUCTION 

NGC 1333 IRAS 4B is a Class young stellar object in 
the Perseus molecular cloud at a distance of 235 pc from 
the Sun (Sandell et al. 1991; Hirota et al. 2008; Enoch et al. 
2009). Radio and infrared observations revealed various star 
forming activities in this region, such as molecular outflows 
and H 2 masers (Blake et al. 1995; Hodapp & Ladd 1995; 
Choi et al. 1999, 2004, 2006; Rodriguez et al. 1999, 2002; 
Choi 2001; Reipurth et al. 2002; Park & Choi 2007). Inter- 
ferometric observations revealed that there are two compact 
continuum sources, BI and BII, with an angular separation 
of 1 1" (Looney et al. 2000; Choi 2001). BI is a protostellar 
object, and most of the star forming activities in this region 
can be attributed to BI. By contrast, BII has been detected 
in the millimeter and submillimeter wavelength bands only, 
and its nature is less clear. 

Lay et al. (1995) performed a single-baseline submil- 
limeter interferometry with a resolution of ^0'/5 and sug- 
gested that IRAS 4BI is a multiple system. However, later 
studies could not confirm the multiplicity of BI (Looney 
et al. 2000; Reipurth et al. 2002). To settle this issue, it is 
necessary to make sensitive observations with a resolution 
higher than 0'.'5. 

While IRAS 4BII has received relatively less attention, 
it is probably an interesting object in understanding low- 
mass star formation processes, because its nondetection in 
infrared suggests that it is an unusually cold object. The lack 
of detectable star forming activity is also intriguing. 

In this paper, we present the results of our observations 
of the NGC 1333 IRAS 4B region in the 6.9 mm and 1.3 cm 
continuum with the Very Large Array (VLA) of the National 
Radio Astronomy Observatory. We describe our radio con- 



tinuum observations in Section 2. In Section 3, we report the 
results of the continuum imaging. In Section 4, we discuss 
the star-forming activities in the IRAS 4B region. 

2. OBSERVATIONS AND DATA 

2.1 Q-band Observations 

The NGC 1333 IRAS 4B region was observed using 
VLA in the standard Q-band continuum mode (43.3 GHz 
or A = 6.9 mm). Twenty-five antennas were used in the C- 
array configuration on 2004 March 2. The phase-tracking 
center was a 20 oo = 03 h 29 m 12?988, <5 20 oo = 31°13'08'/10, 
which is -2" east and -1" north of IRAS 4BII. 

The phase was determined by observing the nearby 
quasar 0336+323 (PKS 0333+321). The flux was calibrated 
by setting the flux density of the quasar 0713+438 (QSO 
B07 10+439) to 0.20 Jy, which is the flux density measured 
within a day of our observations (VLA Calibrator Flux Den- 
sity Database*). Comparison of the amplitude gave a flux 
density of 1.36 Jy for 0336+323. To avoid the degradation 
of sensitivity owing to pointing errors, pointing was refer- 
enced by observing the calibrators at the X band (A = 3.6 
cm). This referenced pointing was performed about once an 
hour and just before observing the flux calibrator. 

Maps were made using a CLEAN algorithm. With a 
natural weighting, the 6.9 mm continuum data produced a 
synthesized beam of ~(X'5 in full-width at half-maximum 
(FWHM). 

2.2 K-band Data 

The NGC 1333 IRAS 4 region was observed using VLA 
in the A = 1.3 cm continuum. Details of the observations 
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Fig. 1. — Map of the A = 6.9 mm continuum emission toward the NGC 1333 IRAS 4B region, made with a natural weighting. The contour 
levels are 1, 2, 4, 8, and 16 x 0.07 mJy beam -1 . Dashed contours are for negative levels. The rms noise is 0.024 mJy beam -1 . Shown in 
the top left-hand corner is the synthesized beam: FWHM = 0753 x 0751 with PA = -1°. The coordinates are position offsets relative to the 
phase tracking center of the Q-band observations. 



Table 1. 

NGC 1333 IRAS 4B Continuum Source Parameters 





Peak Position" 


6.9 mm Flux Density 6 


1.3 cm Flux Density 6 


Source 


a.T2000.0 ^.T2000.0 


Peak Total 


Peak Total 


BI 


03 29 12.009 31 13 08.14 


1.77 ±0.03 3.08 ±0.10 


0.362 ± 0.019 0.47 ± 0.05 


BII 


03 29 12.835 31 13 06.98 


0.77 ± 0.02 0.89 ± 0.05 





a Source positions from the robust- weight 6.9 mm map (Figure 2). The 1.3 cm position of BI agrees with 
the 6.9 mm position within 0'/09. Units of right ascension are hours, minutes, and seconds, and units of 
declination are degrees, arcminutes, and arcseconds. 

6 Flux densities are in mJy beam -1 or mJy, corrected for the primary beam response. 



and the results for the IRAS 4A region were presented in 
Choi et al. (2010, 2011). The observations were made in 
the standard K-band continuum mode (22.5 GHz) in the B- 
array configuration. With a natural weighting, the 1.3 cm 
continuum data produced a synthesized beam of FWHM « 
074. 

3. RESULTS 

Figure 1 shows the 6.9 mm continuum map toward the 
NGC 1333 IRAS 4B region. IRAS 4BI and BII were clearly 
detected. Table 1 lists the continuum source parameters. 
Figure 2 shows a map of the BI region with a smaller beam. 
Figure 3 shows the same region in the 1.3 cm continuum. 
IRAS 4BII was not detected in the 1 .3 cm continuum. 

In the 6.9 mm map (Figure 2), IRAS 4BI seems to be 
elongated in the northwest-southeast direction. An ellipti- 
cal Gaussian fit gives a deconvolved source size of FWHM 



= 0756 x 0'.'15 with a position angle (PA) of -73°. Com- 
parison with the distribution of the H2O maser sources 
(Rodriguez et al. 2002; Park & Choi 2007) suggests that the 
6.9 mm continuum source structure is affected by the out- 
flow, at least in the northwestern part of the source. Possible 
explanations include shock-heated dusts around the outflow 
and free-free emission from ionized outflowing gas. 

In the 1.3 cm map (Figure 3), IRAS 4BI seems to be 
elongated in the northeast-southwest direction. An elliptical 
Gaussian fit gives FWHM = 0720 x 0710 with PA = 50°. 
The size corresponds to 47 x 24 AU 2 at a distance of 235 
pc. Since the elongation is nearly perpendicular to the direc- 
tion of the H2O outflow (PA w 151°), the 1.3 cm continuum 
may be tracing either a circumstellar disk or a flattened pro- 
tostellar envelope (pseudodisk), or probably a combination 
of them. 

The 6.9 mm continuum source of IRAS 4BII (Figure 
1) is essentially unresolved. BII is very compact and has 
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Fig. 2. — Map of the A = 6.9 mm continuum emission toward the 
IRAS 4BI region, made with a robust weighting. The contour lev- 
els are 1, 2, 4, 8, and 16 x 0.07 mjy beam -1 . The rms noise is 
0.024 mjy beam -1 . Shown in the top left-hand corner is the syn- 
thesized beam: FWHM = 0745 x 0740 with PA = -1°. Solid dots: 
H 2 maser sources (Rodriguez et al. 2002; Park & Choi 2007). 
Plus sign: the 3.6 cm continuum source (Reipurth et al. 2002). 
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Fig. 3. — Map of the A = 1.3 cm continuum toward the IRAS 4BI 
region. The contour levels are 1, 2, and 4 x 0.05 mjy beam -1 . The 
rms noise is 0.016 mjy beam - 1 . Shown in the top left-hand corner 
is the synthesized beam: FWHM = 0740 x Of.' 33 with PA = 64°. 
The coordinates are position offsets relative to the phase tracking 
center of the K-band observations. 



never been clearly resolved in interferometric images (e.g., 
Looney et al. 2000; note that IRAS 4BII was referred to as 
IRAS 4C by Looney et al.). 

The total flux densities from our data (Table 1) are 



smaller than those given by Di Francesco et al. (2001) by a 
factor of ~2.5. One possible reason for the difference is that 
our observations were made in more extended array config- 
urations. Another reason could be that the flux calibrator in 
our observations is different from that used by Di Francesco 
et al. (2001). 

4. DISCUSSION 
4.1 IRAS 4BI 

4.1.1 Source Structure 

The single-baseline submillimeter interferometric obser- 
vations of Lay et al. (1995) suggested that the source struc- 
ture of IRAS 4BI cannot be explained as a single source. 
Even a binary model could not fit the data well, and they 
suggested that BI is a triple or higher-order multiple sys- 
tem. Considering their resolution and the clear dips in the 
visibility curve, BI was expected to consist of multiple com- 
ponents of comparable brightness with a separation larger 
than 075. In later studies with resolutions higher than 1", 
however, the image of IRAS 4BI was always found to be 
dominated by a single peak (Looney et al. 2000; Reipurth et 
al. 2002; Choi et al. 2007; J0rgensen et al. 2007; J0rgensen 
& vanDishoeck 2010). 

Our 6.9 mm continuum map (Figure 2) shows a weak 
peak located about 1" northwest of IRS 4BI. This peak is 
probably not a real source, considering that the map shows 
some negative peaks at comparable (absolute) intensities 
and that the 2.7 mm map of Looney et al. (2000) does 
not show a corresponding peak. Even if it is real, it could 
not have caused the interference pattern seen by Lay et al. 
(1995) because it is weaker than the BI main peak by a fac- 
tor of ~9. The elongated shape of the source (Section 3) 
could have caused some interference pattern. Another pos- 
sible source of confusion is BII, because it is located in the 
primary beam of Lay et al. and its existence was not known 
at the time of their analysis. 

4.1.2 Continuum Spectrum and Mass 

Figure 4 shows the spectral energy distribution (SED) of 
IRAS 4BI. (For the SED and mass estimates, we use the flux 
densities given by Di Francesco et al. (2001), 7.8 mjy at 6.9 
mm and 1.2 mjy at 1.3 cm, because the values in Table 1 
might suffer from missing flux.) The steep slope in the (1.3 
cm, 2.7 mm) section suggests that the millimeter flux comes 
from thermal emission of dust. The SED in this section can 
be described using a power-law form, F a where F 
is the flux density, v is the frequency, and a is the spectral 
index. The best-fit power-law spectrum gives a = 3.3 ± 0.2. 
The 3.5-6.2 cm flux may be mostly from free-free emission 
of a thermal radio jet. The free-free emission seems to be 
optically thick at 6.2 cm. 

The mass of the circumstellar molecular gas can be de- 
rived from the dust continuum in the millimeter bands, using 
the procedure and equations described in Section 5 of Choi 
(2001) and the mass emissivity given by Beckwith & Sar- 
gent (1991). The best-fit spectrum gives an opacity index of 
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Fig. 4. — Spectral energy distributions of the IRAS 4B sources. 
Filled circles: BI. Open circles: BII. Flux densities are from 
Rodriguez et al. (1999), Reipurth et al. (2002), Di Francesco et 
al. (2001), Choi et al. (2001), and Looney et al. (2000). The flux 
uncertainties are smaller than the size of markers except for the 
1.3 cm and 6.9 mm data. Dotted lines: Best-fit power-law spectra 
in the (1.3 cm, 2.7 mm) section. 
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Fig. 5.— Maps of the HCN J = 1 -> line toward the IRAS 4B 
region, from the data presented by Choi (2001). For the line core 
of the F = 2 — i> 1 hyperfine component (black contours), the HCN 
line was averaged over the Vlsh interval of 5.8-7.2 km s _1 . The 
lowest contour level and the contour interval are 150 mjy beam -1 , 
and the rms noise is 50 mjy beam" 1 . For the blueshifted and red- 
shifted outflows (blue and red contours, respectively), see Figures 
1 and 2 of Choi (2001) for the velocity intervals. The lowest con- 
tour level and the contour interval are 60 mjy beam - 1 , and the rms 
noise is 20 mJy beam - 1 . Shown at the bottom left-hand corner is 
the synthesized beam: FWHM = 4'.' 6 x 3'.' 3 and PA = 0° . Crosses: 
the 6.9 mm continuum sources. Squares: The H2 line peaks HL 
9a/b (Choi et al. 2006; Figure 6a). 



(3 « 1.3. Assuming a distance of 235 pc to the source and 
a dust temperature of 33 K (Hirota et al. 2008; Jennings et 
al. 1987), the mass of IRAS 4BI is 1.0 ± 0.5 M , which 
may include the inner protostellar envelope and the disk. 
For comparison, the mass estimates of BI in previous works 
(scaled to the distance of 235 pc) are 0.43 M Q in Looney 
et al. (2000), 1.1 Af© in Choi (2001), and 0.47/0.26 M Q in 
J0rgensen et al. (2007). The mass of the whole protostel- 
lar envelope of BI derived from single-dish observations is 
~3.2 M (Enoch et al. 2009, scaled to 235 pc). 

4.1.3 Outflow and Infrared Images 

IRAS 4BI drives a bipolar outflow roughly in the north- 
south direction. Interferometric images of several molec- 
ular lines show the redshifted lobe to the north and the 
blueshifted lobe to the south of BI (Figure 5; Choi 2001; 
Di Francesco et al. 2001; J0rgensen et al. 2007; Choi et al. 
2011). In these millimeter lines the northern and southern 
lobes are often comparable in brightness. Images of the in- 
frared molecular hydrogen line show a jet-like flow to the 
south only (Hodapp & Ladd 1995; Choi et al. 2006). The 
bright H2 knots are referred to as HL 9a/b (Figure 6a). The 
H2 jet may be bipolar in nature, but the unseen northern jet 
may be obscured by the dense protostellar envelope. Recent 
far-IR observations showed highly excited H 2 emission 
from the southern outflow (Herczeg et al. 2011, submitted 
to A&A). 

The direction of the outflow, i.e., the position angle of 
HL 9 with respect to the 6.9 mm continuum source of IRAS 



4BI, is ^174°. By contrast, the position angle of the small- 
scale jet traced by the H2(Dmaseris ^151° (Rodriguez et al. 
2002; Park & Choi 2007). To explain the difference in the 
outflow directions, it was proposed that the jet is precessing 
at a high rate (Marvel et al. 2008; Desmurs et al. 2009). 
(Marvel et al. (2008) also suggested that there are probably 
two outflows, i.e., BI may be a binary system.) The position 
angles of HL 9a and 9b with respect to BI agree within ~ 1 
and show no sign of precession to the proposed direction. 

In the infrared images of the IRAS 4B region, the H2 jet, 
especially the bright knot HL 9a, is the dominating emission 
feature. In the Spitzer images, only HL 9a/b can be seen in 
the 3.6-24 /im bands (Figure 6b-e), and neither BI nor BII 
is clearly detectable. Even at the 70 /im band, the peak po- 
sition is located between BI and HL 9a (Figure 6f), and the 
interpretation of the image and the flux is ambiguous. Note 
that the Spitzer source SSTc2dJ032912.04+31 1301.5 corre- 
sponds to HL 9a, and the SED of this source (Per-emb 13) 
presented by Enoch et al. (2009) is a combination of the flux 
densities of HL 9a on the short- wavelength side and those of 
BI on the long-wavelength side. Therefore, their bolomet- 
ric temperature and luminosity of the BI protostar may be 
overestimated. This example suggests that, at least for some 
young stellar objects, the emission from secondary sources 
(such as outflows and/or binary companion) can "contam- 
inate" the infrared flux and make the SED appear hotter. 
Such effects can introduce a bias into SED surveys of young 
stellar objects. 
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[See http://minho.kasi.re.kr/Publications.html for the original figure.] 

Fig. 6. — IR images of the IRAS 4B region, (a) Map of the H 2 1-0 S(l) line. The emission peak positions are a 20 oo = 03 h 29 m 12? 05, <5 20 oo 
= 31°13'01'.'8 for HL 9a and 03 h 29 m 12!09, 31°12'58'/4 for HL 9b. (b) 4.5 /an image, (c) 5.8 /an image, (d) 8.0 /im image, (e) 24 /an 
image, (f) 70 /mi image. The lowest contour level and the contour interval are 10% of the maximum intensity in each image. The highest 
contour is at the 90% level. The grayscale images are in a logarithmic scale and show low-level features. The H 2 image is from the Subaru 
telescope (see Figure 2 of Choi et al. 2006), and the others are from the Spitzer data archive. 



4.2 IRAS 4BII tinuum (Looney et al. 2000; Choi 2001; Di Francesco et 

al. 2001; J0rgensen et al. 2007; J0rgensen & van Dishoeck 
IRAS 4BII has been detected as a compact source in in- 2010) and in a si le _ dish sub millimeter-continuum map 
terferometric images of millimeter and submillimeter con- 
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(Smith et al. 2000). However, its evolutionary status and 
star-formation activity are not well known. Single-dish ra- 
dio observations usually do not have angular resolutions 
high enough to separate BII from BI, and BII has never been 
detected in infrared. In molecular lines, BII was detected in 
the 13 CO and HCN line maps only (Figure 5; Choi 2001). 
There is no known molecular outflow associated with BII. 
The nondetection in the centimeter continuum (Rodriguez 
et al. 1999; Reipurth et al. 2002) also suggests that there 
is no detectable radio jet. Since there is no detection in in- 
frared, the SED of BII is incomplete and cannot be classified 
well enough to tell the evolutionary status. 

The best-fit millimeter SED of IRAS 4BII (Figure 4) 
gives a mass of ~0.4 M , using the same procedure and 
assumptions in Section 4.1.2. The mass estimate, however, 
is very uncertain because the 6.9 mm flux density (hence 
the opacity index) has a relatively large uncertainty. The un- 
certainty interval of the mass estimate is 0.1-2.3 M Q . For 
comparison, the mass estimates of BII in previous works 
(scaled to the distance of 235 pc) are 0.12 Mq in Looney 
et al. (2000), 2.2 M Q in Choi (2001), and 0.17/0.074 M Q 
in J0rgensen et al. (2007). Note that the dust temperature 
of 33 K used in the mass estimate came from the SED of 
the whole IRAS 4 system (Jennings et al. 1987). Since BII 
is undetectable in infrared, it can be much colder, and the 
mass can be larger than the values given above. 

On the nature of IRAS 4BII, previous studies suggested 
that it may be an optical/IR source (such as T Tauri star) 
or a Class I young stellar object (Looney et al. 2000; Choi 
2001). The nondetection of BII in the Spitzer images (Fig- 
ure 6), however, seems to rule out these possibilities. One 
possible explanation (that was not known at the time of the 
previous studies) is that it can be a newly recognized type 
of young stellar objects known as very low luminosity ob- 
jects (VeLLOs; Young et al. 2004; Crapsi et al. 2005; Lee 
2007; Dunham et al. 2008). The internal luminosity of a 
VeLLO is less than 0.1 L , which is much lower than what 
is predicted by models of low-mass star formation (Shu et 
al. 1987). 

Dunham et al. (2008) found that the 70 ^m flux corre- 
lates well with the internal luminosity. Assuming that the 
strong source in the 70 /xm image (Figure 6f) is IRAS 4BI, 
the image was fitted with a point source. The residual map 
gives an upper limit of ^500 mJy on the flux density of BII. 
Using the scaling given by Dunham et al. (2008), the upper 
limit on the internal luminosity of BII is ^0.08 L Q . There- 
fore, it is possible that BII may be a VeLLO. 

The compact structure of IRAS 4BII and the nondetec- 
tion of any star-formation activity point to yet another possi- 
bility: the first hydrostatic core (FHSC; Boss & Yorke 1995; 
Masunaga et al. 1998). Since the observational characteris- 
tics of FHSC have not been well established yet (e.g., Enoch 
et al. 2010), it is not clear if BII is such an object. More sen- 
sitive and higher-resolution observations in the far-IR band 
are needed to verify the status of BII as VeLLO, FHSC, or 
any other possibilities. 
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